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Abstract: This paper examined the problem of thermosolutal effects in a horizontal cylindrical channel 
with suction velocity. The fluid flow is assumed to be of constant viscosity, thermal conductivity, 
diffusivity and it is axi-symmetrical. The problem modelled followed the work of Ize et al.[24] by 
incorporating suction velocity to their work. Following their method of solution, the governing equations 
were solved. The analytical results were graphically computed using python programming codes and the 
graphs showed that fluid velocity increases with increase heat sink, thermal Grashof and solutal Grashof 
numbers. Also, velocity, energy and species profiles increase with increase in suction velocity while 
thermal conductivity and diffusivity decrease energy and specie profiles respectively. 


Keywords: Thermosolutal, heat, solute, suction, velocity, specie. 


1. Introduction 


Flow in cylindrical channels with suction velocity is of great relevance in natural processes, science, 
engineering, medicine, chemical environments to mention a few. These flow types are influenced by 
coupled heat and mass transport and are applicable in different areas such as oil and gas industries, 
aerospace, chemical and biological processes, heat exchangers, geothermal power plants and groundwater 
movements [1]. Due to these applications, several studies on this flow regime have been carried out by so 
many researchers. For example Isreal-Cookey et al.[2], Nwaigwe [3], Deng and Martinez [4], Rundora 
and Makinde [5], Chinyoka and Makinde[6], Nwaigwe and Makinde [7], Weli and Nwaigwe [8], 
Bunonyo and Amos[9], Emeka et al.[10,11]. 


The combination of energy and specie transport is referred as Thermosolutal transport. Its effect in 
cylindrical channels with suction velocity still have great deal of studies due their relevance in chemical, 
industrial, biological, physical processes and so on. Therefore impact of the thermosolutal transport in 
porous media have been carried out by some researchers such as Hadidi et al.[12] ,Shah et al.[13], Corson 
and Pritchard [14], Mahdy and Ahmed [16], Khan e¢ al. [17], Kefayati and Tang[18]. Nwaigwe[19] 
Futhermore, Nwaigwe [21] numerically solved a heat and mass transport problem in a variable viscous 
channel flow incorporating solute injection. 


In this paper, we extended the work of Ize et al.[23] incorporating suction velocity, Frobenius method of 
power series solution was employed to solve the velocity, energy and specie equations. Our analytical 
results were graphically displayed to capture the effects of the varying parameters on the governing 
equations. 
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This paper is divided into sections as follow: the governing equations including the non-dimensional 
parameters in section 2, while section 3 displayed the method of solution using power series solution of 
Frobenius type in details. In section 4 the graphical results of the analytical solutions are displayed to 
show the effects of the parameters on the governing equations. While the graphical results are discussed 
in section 5. The paper ends in section 6 with conclusion capturing recommendation of further researches. 


2. Mathematical Formulation 


Consider the heat and mass transport in a horizontal cylindrical channel. The flow is considered to be in 
the direction of z° —axis. The fluid’s viscosity and pressure are assumed to be constant. Furthermore, the 
channel walls are stationary with suction velocity. Therefore, the diagram below assumes r’ to be the 
distance of the fluid flow from the centre to the channel wall and v =(0,0,v) is the velocity vector. The 
fluid is also incompressible and Newtonian in nature. 


| 
Stationary Channel Wall 


Tr Ko FLOW AXIS 


Stationary Channel Wall 
Figure 1: The Physical Model and the Coordinate System of the Fluid Flow. 


More conditions imposed are that (i) the fluid flow is of a nonconstant temperature (ii) the flow does 
varies along z’ and it is also steady (iii) the flow is also axi-symmetric. 
Under the assumptions above, the fluid flow is governed by the continuity, energy and specie equations 


(see the previous of Ize et al. [23] for the nondimensional equations) with the following non dimensional 
quantities: 


T -T. C-C. ra Ca ‘ 
= : 2 p= : 2 Gp = SP rlwd Go = Put Ain kt 5 
T, C, UY UY Ko 
* * * 2n* * 2 
peeps ge pet greg Seg i Saye Oe (1) 
a a Uo P% Mo Ce Dy 
T. ° Cpk,T D oP 
rae ) =1- Ne ay ee a aa ~,D= 2 p=-— 
My [, Ko a) D Oz 


Equation (1) are introduced into the dimensional equations to give to the following dimensionless 

equations incorporating suction velocity see Ize et al.[23]. 

Ov _ av ; 1 Ov 
‘ or? rr 


M’v- p+Gr¢-Gco (2) 
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0d Od 10¢ ; 
k = + a 3 
° Or (SS r or G @) 
0p Oo 109 ‘ 
k =D + + 4 
° Or [Es r Or a @) 


Subject to the boundary conditions 


ov _ 06 09 _ 
Or Or Or ° 


(5) 


v=0, &, =%.» %, =P, ONT =1; (6) 


Where &, is the suction velocity, v is fluid velocity component along z’ direction, yz is the viscosity, p is 


the pressure, M is the magnetic field, Gris the thermal Grashof number, Gc is the Solutal Grashof 
number, ¢ is the dimensionless fluid temperature gy is the dimensionless mass concentratio, A is the heat 
sink, @ is the solute injection, « is the fluid thermal conductivity, and D is the diffusivity parameter. 


Since the fluid flow is dependent on z’ only, equations (2) — (6) become: 


dy d’°v ldv 
k,—= + —~M’y— p+Gr¢-G 7 
° dr | 3 r | ae ae ad ) 
do dd 1d¢ 7 
k = + a 8 
° dr {1 r dr é (8) 
dp do ldo . 
k =D + + 9 
° dr [< r ar Oe 0) 
Subject to 
So 0 ae 
or or or (10) 


v=0,¢ =¢,.%, =P, at r=l 
3. Method of Solution 


To solve equations (7)-(10), we assume unique solutions of the forms: 


V(r) = Ayvy(r) + Ay, (7) 


(11) 
PT) = A, (r)+ 4,4(7) (12) 
P(r) = Ay) (r) + 4:9, (7) (13) 


where v,(r), v,(7),@(") . &(7), @(r) andg,(r) are all linearly independent solutions along with 4,, A,, 
A,, A,, A,and A,as arbitrary constants. Therefore, in order to obtain v,(7), ¢(7) and @,(r), Frobenius 
power series soluton is adopted as: ( see Nwaigwe and Amadi [27] and Ize et al.[23] ). 
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Vo (r)=dar"”; (14) 
dy (7) = yr"; (15) 
(16) 


Qo (r) = y &*"; ~~ m=constant 


s=0 
where a,,@&, and d& are all constants to be obtained. 


Therefore, differentiating equations (14) - (16) in their order of degrees in equations (7) — (9) and 
substituting them into equation (7) — (9) with the boundary conditions (10) applied, the results are 


obtained: 
ova (ea [foes +Ma, a} {fies +Ma, a} siamese a} 
(m+2) wu (m+3) (m+4) wu (m+5) yu 
y(r)=r"| 4 k(mt+)a,+Ma, | a4 k(m+Na,+Ma, | , ne k(m+)a+Ma, | oe k(m+)a,+Ma, 0 
: (m+6)" 11 (m+7)" 10 (m+8) w (m+9)° 
| Mant Da, +M 2 ,| Dao +M % |. | holt Daa, +M Ao | 
(m+10)" ye (m+11) Lu (m+12)" u 
(17) 
okie ky(m+) B+ Ak |» oe k(m+)G+A&| 5 ed ky (m+) G4 76 | 4 a Kym Dder 7° B| 
(m+2)" kK (m+3) x (m+4) ( m+5) 
A(ni=r"| + Kym) 2° | | Kolm Der Ake | | mtd 2d | Kolm De 2h | 
(m+6) (m+7) (m+8) (m+9) « 
| Kore Dk 2° 0 | emt +4) | emt +1 &, | , 2 
(m+10)° x (m+ll) « (m+12)’ x 
(18) 


By + e+ K,(m-+1)eRe- of & a | ont Dae Ca oy | hota a | | hole a & . 
(m +2)’ D (m+3) D (m+4)" D (m+5) D 


alr)=r*| +| Seed COME) | Kent Dc AR | Ky (ont ee oC AR) | Ka n+ DoF 
: (m+6)"D (m+7) D 


(m+8) D (m+9)° D 


| here of & Oy k (m+) a — 00 && Ere k(m+1)d&—0f & con 
(m+10) D (m+l1)" D (m+12)’ D 


(19) 
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To obtain v,(r), ¢(r) and g(r) in equations (14) — (16), equations (17) — (19) are differentiated with 


respect to m to obtain: 
pM ecie ee 


wr vael 
k,(m+l)a, +Ma, k,(m+l)a, +Ma, k,(m+l)a, +Ma, 
(m+7) uw (m+8) wu (m+9) uw 
k,(m+)a, +Ma, 


‘ (m+10)" y2 


k,(m+)a,+Ma, 
(m+2)° wu 


re 
H 


k,(m+)a,+Ma, 
(m+3) ys 


re 
Hl 


k,(m+la,+Ma, 
(m+4) w 


rt 


k,(m+la,+Ma, 
(m+5) y0 


y 


| font Na, +Ma, | « 
(m+6) uw 


rior 


k,(m+])a, +Ma, 
2; 
(m+11) L 


k,(m+)a,,+M ay 
(m+12) L 


(p= 


k,(m+la, +M M 


ay” Fart (i 


k,(m+la,+Ma, 


k,(m+l)a,+M* a 


k,(m+l)a,+Ma, 


(m+2) wu 


re 


(m+3) wu 


al 


(m+4) wu 


rt 


} 


(m+5) wu 


k,(m+la,+Ma, k,(m+l)a,+Ma, k,(m+)a, +M'a, | , 4A 


o(m+1)a, +.M? a 


pein FC 


rt 


(m+6) wu (m+7) wu (m+8) wu 


k,(m+l)a, +M Ay 
(m+12) yu 


k,(m+Da, +M ak 
(m+10) 


k,(m+))ay + Ma, 
(m+11) w 


t 
rl 


p 


(m+9) 


(20) 


0 ge | Kom De Ak 24 
4 re (m+2) « | 


+ ba) : | 


ky ad ay S 


ky(m+ B+ Vb 
(m+3) « 


(m+4) « 


r’lnr 


+6) K 


+1)&+ Vk 04 
(m +10)’ x 


k, ed * 


Kk (m+) b+ Vb | 
(m+7) K 


(m+8) K« 
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(m+12) 
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Ja c 
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Jed) 


k(m+ e+ - ba * 


ky (m+ b+ 2.8 |; 3 
(m+5) 


| 
| 


+ Ab a, 
(m+6) K 

Ky (m+ Pe) 
(m+10) x 


k, ed is 
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ky (m+l)& + 1b ay 
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(21) 
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fi’ ve! (n+ 1) a 24); {i (n+ 1) i {: (n+ 1)8—ar 28, {: (n+ 1) 07 a, 
(m +2)'D (m+3) D (m+4)' D (m+5) D 
pr] | lem Dao |, _[ Kolin De 07 | | ky m+ Dd 07) | k(t De-0| 
(m+6) D (m+7) D (m+8) D (m+9) D 
[Seen co | ine 8 | See ele 
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| (m+2)° D (m+3) D (m+4) D (m+5) D 
i é F k(m+]) dar | ky (m+ )&&-a? a 44 (m+ 1)d&— a? a | s 44 (m+ |) ar a) 5 
am (m+6) D (m+7) D (m+8) D (m+9)' D 
{Ral Di 07 AE) | (on + Nd — 7A | ky n+ dR — 0 | 
(m+10) D (m+1l) D (m+12) D 


(22) 


Substituting for m=0 and writing the equations (20) — (22) in terms of a,, we obtain the velocity, 
temperature and concentration solutions as: 


(My (2 } (6M +9M" wr" ,{ CAMs +68M'k pi 
Qu Leye (4!) ue (5!) 4é 
{ones + 400M 1+ 25M p25 } [came +3804 4 378M Kp yy" 
6! 1 “ 
ur)=4, ( ) e ( ) e +a ae v(r))+ ae 


(5040M!°K> +32508M"k;, piesa Le +11025M*k pe yr 


{ 


(81) a 
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; 


{ 


(9a 


I 


2 3 4 5 6 7 8 9 
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(ary (ae (62 +9209" | | (242 +0824 wr? 
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(2a (GBI (4!) 1° (51) x4 


{ms $A4NA FE +DEA WS = $3804.04 37984 Ky" 


6l) TI | 
(6) (7!) A sf 4on)+08 
| (S04027K} + 325082} -+505962°K A? +110052° Ka an 
(8!) x’ 
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A 
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, (ary Gaol (607K; +9at'Dyr* |_| (240°? + 680k, Dyr* 
(2!°D \G!D (4!) D (5!) D* 
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(6!) D’ (7!) DS 
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om ) 
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(9!) D’ ~ 
(25) 


Hence, the unique solutions for velocity, temperature and concentration as expressed in equations (11) — 
(13) are obtained as: 


ay ee (6M +9M ir" \ | 24M +68 K pi? 
a wo) a 2 + 2 + 2 
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(26) 
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(28) 
Putting the boundary conditions in equation (10) into equations (26) — (28) and setting 4, = A, = A, =0 
due to boundedness the following results are expressed as: 
2 2 24,2 4 
ae. A or} «(2h An wn, Je es ea HA, +T, |r? 
(2!) us (3!) u (4!) 
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Ary { 2k,4?r° (617k +9A%K«)r* (2417k° + 681*k,.K)r° 
+ +| —° + : + u : 
(2I°« \GBYK (4!) «° (5!) x? 
(120A*kS + 490A *k +2251’ )r® (720A°k. +380417*k}K+3798A°k Kr’ 
Se OS pe ee 
2 2 
ee d, (6!) « (7!) «° 
#(D) _[ (S040A7k) +32508A*kjx+50596A%K5K* + 110252 kya)" 
(8!) x” 
(40320A7k) +3061441*koK + 6482241 °k3K? +331272A5k, «er? 
+ 0 o 0 0 2 oe 
(9!) «* 
(30) 
= (ary _{ 2k,arr’ 7 (-6a°k; + 9a*D)r* : (—24a7k; + 68a*k,D)r° 
(2!)°D ((!)°D? (4!) D® (5!) D 
_ (-120a7k} + 490a*k; D-225a°D’)r® |_| (-720a7k; +3804a*k; D-3798a°k,D’ yr’ 
acer er) (6!) D° (7!) D’ 
7 (-504007k§ +32508a*k; D —50596a°k; D? +11025a°k,D’)r* P 
(8!) D’ ~ 
(31) 
where A, = e, , and A, = Pw 
%,, (1) 9,, (I) 
4. Results 


The graphical results of the analytical solutions for flow velocity, temperature and sepcie solutions are 
displayed in this section. The three solutions are plotted against different parameters such as 
viscosity,thermal Grashof, solutal Grashof, magnetic field, thermal conductivity, heat sink, solute 
injection and diffusivity at varying values, so as to show the effects of suction velocity. 


velocity 


0.0 0.2 0.4 0.6 0.8 1.0 


Figure 2: Velocity Profile against r for the values of Gc=0.5, Gr=0.5,a=0.6,k=0.3, D=5,M=5, kp =10, 
A=2 varying viscosity parameter. 
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Figure 3: Velocity Profile against r for the values of Gc=0.5, Gr=0.5, o=0.6, k=0.3,D=5, M=5, kp=0.3, 
p=0.1, p=1.0, varying heat sink parameter 


velocity 
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Figure 4: Velocity Profile against r for the values of a=0.6,K=0.3.D=5,M=5, Gr=0.5, p=0.1,A=2, 
ko=0.3,p=1.0 varying Thermal Grashof Number 


velocity 


Figure 5: Velocity Profile against r for the values of Gr=15,a=0.6,k=0.3,D=5,M=0.6, kp=0.1, »=0.1, 
p=1.0, varying Solutal Grashof Number. 
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Figure 6: Velocity Profile against r for the values of Gce=0.5, Gr=0.5,K=0.3,D=5,M=5, p=0.1, 
A=0.8,p=1.0, varying suction parameter. 
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Figure 7: Temperature Profile against r for the values of 0w = 0.1, 4 = 0.1, ko= 0.1, varying thermal 
conductivity parameter. 


Temperature 


Figure 8: Temperature Profile against r for the values of 0w = 0.1, k = 1.0, ko= 0.1, varying heat 
sink parameter. 
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Temperature 


Figure 9: Temperature Profile against r for the values of Ow = 0.1, kK = 1.0, 4 = 0.1, ky = 0.1, varying 
suction parameter. 
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Figure 10: Concentration Profile against r for the values of Cw = 0.1, a = 1.0, ky=0.1, varying 
Diffusive parameter. 
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Figure 11: Concentration Profile against r for the values of Cw = 0.1, D= 1.0,K0=0.1, varying solute 
parameter. 
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Figure 12: Concentration Profile against r for the values of Cw = 0.1, a = 0.1,D=5.0, varying suction 
parameter. 


5. Discussion 


The problem of thermosolutal effects have been solved extending our previous work in Ize et al. [23] with 
the inclusion of suction velocity and it is discovered that, to a large extent, the inclusion has added more 
value to the work. Therefore, the discussions coming from the study are based on heat sink, solute 
injection, diffusivity, thermal conductivity, thermal Grashof, Solute Grashof numbers and suction velocity 
which are our main interest. The discussions are done with the help of graphs and the varied values of the 
parameters in the problem are indicated in the graphs. 


Figure 2 shows the effect of viscosity on the fluid velocity. It is observed that flow velocity increases 
steadily close to the wall as viscous parameter increases. This is due the fact that suction reduces frictional 
losses in fluid thereby enhancing fluid velocity see Kishan and Amrutha [15]. In Figure 3, the flow 
velocity shows no significant change at the centre of the channel as the heat sink is varied but as it moves 
away from the centre to the wall, velocity increases with heat sink increases. This is in agreement with 
Asogwa et al. [22] 


In Figure 4 and 5 the effects of thermal and solutal Grashof numbers on the velocity profile of the flow 
are shown and observations show that they both increase the fluid velocity as they are varied see 
Nwaigwe [19]. 


Figure 6 shows the effect of suction parameter on the fluid velocity distribution. It is discovered that a 
surge in suction parameter increases the flow velocity and this is more pronounced at the wall where 
suction reduces the drag of the body of the flow. This is in agreement with Mebine & Adigio [20]. The 
effect of thermal conductivity parameter on temperature profile is illustrated in Figure 7 and it is 
discovered that fluid temperature decreases as the thermal conductivity parameter increases. Also Heat 
sink parameters have effect on fluid temperature in figure 8. It is realised that the fluid temperature 
increases more at the wall as the cooling parameter increases but no effect is seen at the centre. In Figure 
9, the influence of suction on fluid temperature is shown. It is observed that fluid temperature increases at 
the wall as suction parameter increases. 


The significance of diffusive parameter on fluid concentration in figure 10 is displayed. It is discovered 
that the fluid concentration decreases as the diffusive parameter increases. While Figures 11 illustrates the 
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effect of solute injection parameter on fluid concentration. Observably, increase in the solute injection 
term depicts increase in the fluid concentration. Also, the suction parameter has an impact on the fluid 
concentration in figure 12. It is discovered also that, as the suction parameter increases the fluid 
concentration. 


6. Conclusion 


In this study, thermosolutal effects in cylindrical channel have been studied and the analytical solutions 
for velocity, energy and specie equations are displayed graphically using python programming codes. 
Results showed that fluid velocity increases with increase in parameters such as heat sink, thermal 
Grashof and solutal Grashof numbers. Also, suction velocity increases the velocity, energy and species 
profiles while both thermal conductivity and diffusivity decrease energy and specie profiles respectively. 
Further studies can be done by incorporating variable viscosity and thermal conductivity see Nwaigwe 
[24]. In this case a stable and convergent numerical method may be needed in both one dimension see 
Nwaigwe [24,26] and two dimensions [25]. 
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